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Abstract 

By resorting to sequential constructions of exchangeable random partitions (Pitman, 2006), 
and exploiting some known facts about generalized Stirling numbers, we derive a generalized 
Chinese restaurant process construction of exchangeable Gibbs partitions of type a (Gnedin 
and Pitman, 2006). Our construction represents the natural theoretical probabilistic framework 
in which to embed some recent results about a Bayesian nonparametric treatment of estimation 
problems arising in genetic experiment under Gibbs, species sampling, models priors. 
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1 Introduction 

In two recent papers (Lijoi et al. 2007, 2008) a Bayesian prior to posterior analysis for the subclass 
of exchangeable partitions in Gibbs form of type a € (0, 1), first introduced in Pitman (2003), and 
largely studied in Gnedin and Pitman (2006), has been proposed for a nonparametric treatment of 
some inferential problems arising in genetic experiments. 

Here, we derive a generalized group sequential construction of exchangeable partitions in Gibbs 
form of type a to place this theory in its natural probabilistic framework and to provide new 
insights on the derivation of relevant results. Our construction, which relies on known results of 
Pitman (2003, 2006) and Gnedin and Pitman (2006), has potential applications for investigating 
additional distributional results for quantities of statistical interest when exploiting in a Bayesian 
nonparametric perspective the theory of exchangeable partitions. 

Notice that, while in Lijoi et al. (2007, 2008) the treatment is in terms of factorial coefficients, 
our treatment is in terms of Stirling numbers, which naturally arise when summing over spaces of 
partitions with fixed number of blocks (see e.g. Pitman, 2006, Ch. 1) and also admit a generalized 
version. Moreover many convolution, multiplicative and multinomial formulas hold for Stirling 
numbers that greatly improved and simplify the presentation. 

The paper is organized as follows. In Section 2 we recall some preliminaries and basic definitions 
on rising factorials, random partitions and generalized Stirling numbers. In Section 3 we focus on 
infinite exchangeable partitions and derive a group sequential version of the Chinese restaurant 
process for Gibbs partitions of type a. In Section 4 we show how to embed some results in Lijoi et 
al. (2007, 2008) in our setting. 
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2 Preliminaries and basic definitions 

We start by recalling some known facts about rising factorials and Stirling numbers which we will 
largely exploit in the following. A comprehensive reference for the role of these numbers in the 
theory of exchangeable random partitions is Pitman (2006). For the sake of clarity we strictly 
adopt his notations. 

For n — 0, 1, 2, ... , and arbitrary real x and h, let (x) n -\h denote the nth factorial power of x 
with increment h (also called generalized rising factorial) 

n-i 

{x) n]h := x(x + h)---(x + (n- l)h) = JJ (x + ih) = h n (x/h) nh (1) 

i=0 

where (x)„| stands for (x) n ^i, and (x)o-\h = x , for which the following multiplicative law holds 

{x) n+ r}h = {x) n]h (x + nh) r 1h- (2) 

From e.g. Normand (2004, cfr. eq. 2.41 and 2.45) a binomial formula also holds, namely 

{x + y) nVl = [ 7 !)(x)^h(y),i-kTh, (3) 
fc=0 ^ " ' 

as well as a generalized version of the multinomial theorem, i.e. 

v , v 

E*J'W»= „>'['.„ ; \ i ( 4 ) 

j=l nj>0,J2 n j= n ' P ' j = l 

Notice that for rrij > 0, for every j, and Y]j rrij — m, an application of ^ yields 

{Zj) nj+mj -l1 = {Zj) m ^ h {Zj + TUj - l)„ iT (5) 

and by (jH) 

, p p p 

/fcl . ' ' ' / bir\ • 

n 3 >G,Y.nj=n P j=l j=l j=l 

P P 

3 = 1 j=l 

which agrees with the result in Lemma 1. in Lijoi et al. (2008). 

Now recall that a partition of the finite set [n] — (l,...,n) into k blocks is an unordered col- 
lection of non-empty disjoint sets {A, . . . , A^} whose union is [n], where the blocks Ai are assumed 
to be listed in order of appearance, i.e. in the order of their least elements. Recall also that the se- 
quence (| Ai|, • • • , | At |) of the sizes of blocks, (m, . . . , n^), defines a composition of n, i.e. a sequence 
of positive integers with sum n and call VfL the space of all partitions of [n] with k blocks. 

From Pitman (2006, cfr. eq. (1.7)) we know that the number of ways to partition [n] into k 
blocks and assign each block a W combinatorial structure is given by 

fc 

B n<k (w,):= Y Il^-I ( 6 ) 
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which is a polynomial in variables w, — (wi, . . . , w n -k+i), known as the (n, k)th partial Bell polyno- 
mial. Now, for each unordered partition of [n] into fc disjoint non empty blocks there are fc! different 
ordered partitions of [n] into k such blocks, and corresponding to each composition {n\, . . . , nk) of n 
with k parts, there are ( ™ n ) = n\ Yli—i different ordered partitions {A\, . . . , -Afc) of [n] with 
= n,-. Hence the definition of B n> k{w») as a sum of products over "Pt^, translates in terms of 
sum over compositions of n into fc parts as 

(ni,...,7ifc) «=1 

In what follows we will largely exploit the notion of generalized Stirling numbers, (for a compre- 
hensive treatment see Hsu and Shiue, f998; see also Pitman, 2006, eq. 1.19). For arbitrary distinct 
reals 77 and (3, these are the connection coefficients S^'f defined by 

n 

( X )n.lr, =Y. S n.k( X )klP 

and correspond to 

For 77 = —1, f3 — —a, and a G (—00, 1), S~ is denned by 

n 

0)«T1 = E S njc~ a ( x )k1 a - (7) 

For w ni = (1 — a) ni _i| and a € [0, 1), equation ([6]) yields 

im(i-«)._ 1T )= e na e n ' 1 '; 1 :- 1 ^ ^ w 

{Ai,...,Aj,}eP [n] k*=l (ni,...,n t ) i=l 

Remark 1. As recalled in the Introduction, in Lijoi et al. (2007, 2008) the treatment is in term of 
generalized factorial coefficients, which are the connection coefficients C(n, fc; a) defined by 

n 
k=0 

(cfr. Charalambidcs, 2005). From ([T]) and |(7J), if x = ya then 

n n 

(ya) n n = E S n-k~ a (y a )kta = E ^aT""^)^ 
k=0 k=0 

hence 



not 

The representation (37) in Lijoi et al. (2008), (Toscano, 1939), also holds for generalized Stirling 
numbers with the obvious changes (cfr. e.g. Pitman, 2006, eq. 3.19). Additionally, specializing 
formula (16) in Hsu and Shiue (1998), the following convolution relation holds, which defines non- 
central generalized Stirling numbers 



s=k 
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and by ©, 

= a k S^ = (")cj; fc (-7)«-.Ti. 

Hence the following variation of equation (38) in Lijoi et al. (2008) defines non-central generalized 
Stirling numbers as connection coefficients, 

n n 

{ya - 7)n?i = Y S n^ a ' lak iy)m = Y s n^~ a ' 7 (y a )kt»- ( n ) 

fc=0 k=Q 

3 Generalized sequential constructions of exchangeable par- 
titions 

Random partitions are random objects that arise in many contexts, exchangeable random partitions 
arise e.g. by sampling from random, almost surely discrete, probability measures. First recall that 
given a law Q on the space v\ of decreasing sequences of positive numbers with sum 1, and a law 
H(-) on a Polish space (S,S), a random discrete probability measure (RDPM) P on S may always 
be defined as P( ) = Y^iLi Pi&Xi {')■> f° r X-i ~ H(-) and (Pj) ~ Q. From Kingman's theory of 
exchangeable random partitions (Kingman, 1978), sampling from P induces a random partition II 
of the positive integers N by the exchangeable equivalence relation i rj j -O- Xi = Xj , that is to say 
two positive integers i and j belong to the same block of II if and only if Xi = Xj, where Xi\P are 
iid ~ P. It follows that, for each restriction n„ = {A\, . . . , A k } of II to [n] = {1, . . . , n}, and for 
each n = 1, 2, . . ., 

Pr(II„ = {Ai, . . . ,A k }) =p(ni, . . . ,/ifc), 

where, for j — 1,2, . . . , k, nj — \ Aj\ > 1 and 53j=i n j — n i f° r some non-negative symmetric func- 
tion p of finite sequences of positive integers called the exchangeable partition probability function 
(EPPF) determined by II (see Pitman, 2006, for a comprehensive account on exchangeable random 
partitions and related stochastic processes). 

The particularly tractable class of exchangeable random partitions of Gibbs form of type a has 
been first introduced in Pitman (2003) and then studied in Gnedin and Pitman (2006). Results 
and explicit forms for many of its EPPFs have been recently obtained in Ho et al. (2007). Here we 
recall the basic notions. An exchangeable random partition n„ of the first n positive integers, is 
said to be of Gibbs form if for some nonnegative weights W = (Wj) and V = (V n ,k) the EPPF of 
n can be expressed in the product form 

k 

p(n 1 ,...,n k ) = V n ^Y[W nj (12) 
j=i 

for all 1 < k < n, and all compositions (ni, . . . , Uk) of n. Gnedin and Pitman (2006) show that to 
define an infinite random partition of N, i.e. a sequence (n„) consistent as n varies, satisfying 

fc(n)+l 

p(n) = p(ni, . . . , n k ) = Y p(ni,---,nj + l,...,n k ) (13) 

j'=i 

for all compositions of n, the weights (Wj) must be of the form W nj = (1 — a) n ._if for a € [-co, 1), 
(with Wj — 1 for every j, for a — — oo), and the weights (Vn.fc) must be the solution to the backward 
recursion V n ,k = (n — ak)V n +i t k + V^+i^+i with Vi,x = !■ They also obtain the solutions, for each 
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a < 1, identifying the extreme points of the infinite dimensional simplex of the possible weights V, 
and deriving corresponding families of extreme partitions, in terms of the laws of the corresponding 
ranked atoms (Pi). Their fundamental result, already stated without proof in Pitman (2003, cfr. 
Th. 8), is the following: 

Theorem 2. [Gnedin and Pitman, 2006; Th. 12] Each exchangeable Gibbs partition of a fixed 
type a S [— oo, 1), i.e. characterized by an EPPF of the form 

k 

p{ni, ...,n k ) = V n ,k - aOr^-rr 

3=1 

is a unique probability mixture of extreme partitions of this type, which are 

i) for a G [— oo, 0) PD(a, m\a\) partitions with m = 0, 1, . . . , oo, 

ii) for a = PD(0, 9) partitions with 9 G [0, oo), 
Hi) for a £ (0, 1) PK(p a \t) partitions with t £ [0, oo). 

Recall that PD(a, 9) stands for the two-parameter Poisson-Dirichlet distribution (Pitman and Yor, 
1997) and PK(p a \t) for the conditional Poisson-Kingman distribution derived from the stable sub- 
ordinator, (cfr. Pitman, 2003). 

Now it is known (cfr. Pitman, 2006, Ch. 3), that the consistency condition 1|13|) allows to derive a 
sequential construction of exchangeable partitions, known as Chinese restaurant process construc- 
tion (first devised by Dubins and Pitman in 1986), that, in its more general formulation, (cfr. 
Ishwaran and James, 2003) is usually introduced as follows: 

Given an infinite EPPF, p(n) — p[n\ , . . . , nj,), assume that an unlimited number of customers ar- 
rives sequentially in a restaurant with an unlimited number of circular tables, each capable of sitting 
an unlimited number of customers. Let the first customer to arrive be seated at the first table. For 
n> 1, given (m, ... ,nk), the placement of the first n customers at k tables, the n+lth customer is: 

a) seated at the table j, for 1 < j < k n , provided p(n) > 0, with probability 

p(n j+ ) 

Pj,n=Pj(n)= (15) 

where p(n J+ ) stands for p(n\, . . . ,rij + 1, . . . , nk) or is 

b) seated at a new table with probability 

Pa.n = Po.n(n) = — -— (16) 
p(n) 

for I = k n + 1, and Y,jiiPj,n + Po,« = L 

Here we derive a group sequential version of the Chinese restaurant process which can be intro- 
duced as follows: 

Given an infinite EPPF p(n), assume that an unlimited numbers of groups of customers arrive 
sequentially in a restaurant with an unlimited numbers of circular tables, each capable of sitting 
an unlimited numbers of customers. Given the placement of the first group of n in a (ni, . . . , n^) 
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configuration in fc tables, the new group of m customers is: 

a) seated a' 
probability 



a) seated at the old fc tables in configuration (mi, . . . , TO&), for mj > 0, 5Z,- =1 Tnj — m, with 



*»(») = p(m|n) = ± ^ ' ' ' ' n \ + ^ , (17) 

b) seated at fc* new tables in configuration (si, . . . , Sfe*), for X)y=i s j = m, 1 < fc* < m, Sj > 1, 
with probability 

/ \ , I vP(ni,---,nk,3i,...,3 fc .) 

p B (n )=p(bdl) 7- , 18 

p{m, ■ ■ .,n k ) 

c) s < m customers are seated at fc* new tables in configuration {s\, . . . , Sfc*) and the remaining 
to — s customers at the old tables in configuration (toi, . . . , rrik) for X^JLi m j = m — s, 1 < s < m, 

j Sj — s, rrij > 0, Sj > 1 with probability 

Pm , s (n) = P (m, s|n) = y + (19) 

p(ni,...,n fc ) 

From now on we focus on the particular case of Gibbs EPPFs. The mathematical tractability of 
the Gibbs product form, combined with the properties of generalized rising factorials previously 
recalles, allows an easy derivation. 

Proposition 3. For the Gibbs EPPF iflij). formula (IT). (HI) and {JTHI) specialize as follows. 
Given the placement of the first group of n customers in a (m, . . . ,n k ) configuration in k tables, 
the new group of m customers is 

a) seated at the k old tables in configuration (mi, . . . , m^), for mj > 0, Sj=i TO i = m : with 
probability 

/ s V n+m,k ]lj=l( 1 - Ot) n +m l 

Pm{n) = v rrfc n ^ (20) 

^n,fcll J= i(l ~ a)„ 3 _i T 

which, by means of the multiplicative relation [[5]). simplifies to 



r, 



n+m,fc 



b) seated at fc* new tables with configuration (si, . . . , s/c*), for Y]j—-\ Sj — m, 1 < fc* < to, Sj > 1, 
with probability 

at » K+m.fc+fc* Il^ll 1 - a )"i-lT Il^ll 1 ~ «)sj-1T V„+m,fc+fc* ■fr, , . . 

P( n ) = nfc 7j : = 7TT 11C 1 ~ «) Sj -i T , (22) 

c) a subset s < to of the new customers is seated at fc* new tables in configuration (s±, . . . , s^*) 
and the remaining to — s customers are seated at the old tables in configuration (mi, . . . ,m k ) for 

m j = m — s, 1 < s < m, X^=i s j ~ s > m j — 0, Sj > 1 with probability 

N V n+m ,k+k* IlLlC 1 " a)n 3 +mj-lT IlLlC 1 ~ a )s 3 -lT 

Pm( n ) = ^ = 
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which simplifies to 



n-\-m,k-\-k 



1 I[(^-«) mjT n(l-a) S3 .- 1T . (23) 



V n .k 

3 = 1 J=l 



Corollary 4. Consider the event A={A11 m new customers are seated at new tables}, by (JH]) and 
(|22|) . summing over all the way to allocate m in fe* tables for every k* yields 

fe*=l Sj>l,£j s 3= m V 7 ' 3=1 

m 

E yn+m,fc+fc* q— 1,— a 
f7~ m,fc* ■ 

fc* — 1 ' 

For the event B={A11 m new customers are seated at the k old tables}, by (|2"Tj) it is easy to obtain 

^< B »« ->=T?f , e , L™ Jrk- -«>-,.- <*> 

(mi,..,m fc ) v 7 j=l 

£\ rrij— m,mj>0 

and specializing (|4j) for #j — (rij — a), 

= 77 (>>j -a))mT = -T7 (n-fea) roT . 

3=1 

4 Embedding the Bayesian nonparametric analysis in the 
group sequential construction. 

Now we show how to embed in our construction some of the results in Lijoi et al. (2007, 2008). 
Notice that our derivation shows explicitily how generalized Stirling numbers arise in this context, 
hence indirectly clarifies how generalized factorial coefficients arise in their treatment. 



Proposition 5. The result in Proposition 1. of Lijoi et al. (2008) may be obtained by 
by summing over the ways to choose (m — s) integers from m, and all the ways to allocate (m — s) 
integers in the k old tables. 

Proof. Marginalizing (f23|) with respect to (mi, . . . , mk) yields 

p(s|n) =p(si ) ...,s k *\ni,...,n k ) = 

k k* 



V n ,k \m-sj ^ \m 1 ,...,m k/ 



e l, m jn^--wri( i -^T= (26) 

,...,m k ) 

■ m j — rn -s,mj>0 



and by (0} 



Vn+m,k+k* ( m 



V n k \m - S 



(n - ka) m - s ^ rj(l - a) Sj ._i T . (27) 
j'=i 
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Corollary 6. Corollary 1, in Lijoi et al. (2008, eq. (10) and (11)), may be obtained from 

^k* 
J 3~ 



27|) by summing over the space of all partitions of s elements in k* blocks. For s = Ylj=i s ji an 



application of (TSJ) yields 



Pr(K* =k*,S = e\m, ...,n k )= V »>+™> k + k ' f m ) („ _ ka) m ^ . (28) 

*n,k V s / 

Marginalizing with respect to k* a probability distribution for the total number S of observations 
in new blocks is easily derived in terms of generalized Stirling numbers, 

Pr(S = s|ni, . . . ,n k ) = — — ( m J (n - /ca) m _ sT ^ Vn+m,k+k* S~l',~ a . (29) 



fc*=0 



Proposition 7. [Lijoi et. al (2007, eq (4)] The probability distribution of the number of new 
blocks k* may be derived from (f^5|) and exploiting the convolution formula (JTOJ) . 



Proof. By marginalizing (|28p with respect to S, 

Pr(K* = k*\ni, ...,n k )= Vn+r "' k+k ^ f m j („ _ fca) m _ sT S'^r Q . (30) 
From (TTU1) wc know that 



S»,fc~ a ' 7 - E { s ) s s,l'- a (-^n- s n (3i) 

hence 

Pr(JT = fc*|nx, . . . ,n fc ) = Vn +™> k + k ' S m \- a ^ (n - ka) , (32) 

where S' m 1 ^,. ^™ fea ' is the non-central generalized Stirling number defined in pip for scale pa- 
rameter 7 = (n — ka). The result in equation (12) in Lijoi et al. (2008), expressed in terms of 
Stirling numbers, follows from the ratio of (f2"8"|) and 



( m ) (n — ka) m -itS I', 
Pr(S = s\K* = k\ ( ni) . . . , n k )) = ^ >- * T s > h 



q— X,— a, — (n— hoc) 
^m,k' 

By interpreting distributions (|2"9"|) and (|32|) as posterior distributions in a Bayesian perspective the 
corresponding expected values give corresponding Bayes estimator under quadratic loss function 
for K* and S, (see Lijoi et al. 2008, eq. (14) and (15)). 

Conjecture 8. In Proposition 2. Lijoi et al. (2008) obtain the following simplification for the 
expected value of the number S of observations in new blocks, given the placement of the first n 
customers 

E(b\ni, ...,n k )= m — — , 

V n ,k 

and provide a proof by induction based on properties of conditional expectations under exchange- 
ability 



8 



In our setting we conjecture that an alternative proof may be obtained by suitably exploiting the 
backward recursive equation for the Gibbs weights 

V n ,k = (n - ka)V n+ i t k + 
combined with some already recalled relationships for Stirling numbers. 

First notice that by ©, for positive x and y 

m\ (x) s (y) m - s 



^ ' c ' (x + y) m 



s=0 

and by substitution z = x + 1 and t = s — 1 it easy to show that 

fm\ {x) s {y) m - s _ / x \ /m - 1\ (z)t(y) m -i- t 
(x + y) m "V^ + y/^V « / (« + tf)m-i 

Therefore, if one is able to show that the probability distribution of S 1 



/V( ,S' = .s| /( i ;/ A ) = -j-i- ( m j (n - fca) m _ sT 14 +m)fe+ fc» 5 S ^; Q (34) 

fc*=0 



V n ,k V S 

can be reduced to 



Pr(5 = S K, . . . , n k ) = ( m ) (("-MK + i, fc ) m -T(K + i, fc+ i). T (35) 
V s / (K,fc)m T 

the result would follow. We conjecture that the intermediate step would be to prove that the 
backward recursive equation sufficies to show that (|34[) equates 

s 



)(»-fea) ro _ sT n + m < k ^ (V n+ i,k+i)k*TaS St l* (36) 

/ V n,k k , =Q 



in fact, by means of ([TJ) , it would reduce to (|35p . We remark that, up to now, this is a conjecture 
and equivalence of (J34J, (|36f and (|35|) is still to be proved. 

Example 9. The previous conjecture is inspired by the particular case of the two parame- 
ter Poisson-Dirichlct model (Pitman and Yor, 1997). It is known that it belongs to the Gibbs 
class of type a £ (0,1) since arises by mixing the stable conditioned PK(p a \t) model by j(t) = 
r(e/«+i) t~ 6 f<* ft) ■ The Gibbs weights are known to be 

v _ {8 + a)k-ita 
Vn ' k ~ (0+l)„_ 1T ' 

By exploiting the recursive equation some algebra shows that 

K+i,fc+i _ 1 _ (n - fca)K+i,fc _ (0 + fca) 
V n ,k ~ V n , k ~ (0 + n) 

and 

V nt k (0 + n) m ' 
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By © and (J7J), the sum in ||25J) results 



1 (a . •> a-l-a {9 + a)k-l-\a , , a -l,-a {0 + a)k-lla (Q . , x 

— — ^ (fl+a)k+ fc ._i Ta 5 fc : = TT— - t 2^ ( e +Mfen« 5 s .fei = 7^—7^ {0+ka) s1 . 

(O + lj n+m _i (0 + lj„+rn-l (V + IJn+m-l 

hence equation (|29|) specializes as 
Pr(5 = s|ni, . . . ,nk) = ]{n - ka) m _ s ] — — '—(9 + ka) s |— ' ' 



and by (f33|) the expected value results 

£7(6 Tlx, . ..,7ik) = to— — . 

ft + n 

Corollary 10. [Proposition 4, Lijoi et al. (2008)] The probability that the m new customers 
don't seat at a subset of (k — r) old tables arises from (f2"3"| by summing over the ways to choose s 
customers from the m of the new group, by summing over the ways to partition s customers in a 
subset of k* new tables, for k* > 1, and over the ways to allocate m — s customers in at most r old 
tables. Notice that the operations of partitioning and allocating differ for the fact that the blocks 
of partitions cannot be empty while allocation in a fixed number of blocks can result in a certain 
number of blocks remaining empty. From (|23p we obtain 

fc*=l n ' K s=k* V 7 s u ...,s k *,-E. Sj =s V " 7 3=1 



X 



(mi ,. . .,77V) m j = m — s J — 1 



which simplifies to 



K fc EL) (E ^ - ™)rn-s T 

fc*=l ™' K s=fc* V 7 j=l 



and the definition of non central generalized Stirling numbers (cfr. eq. (I10p ) yields 

m T ^ 

, =1 m ' fc * 



4.1 The reproducibility of the Gibbs structure from "deletion of classes" 
property of PK models 

In Section 3.1 Lijoi et al. (2008) point out a "reproducibility" property of the exchangeable Gibbs 
partitions of type a which motivates the definition of conditional Gibbs structures. Here we investi- 
gate the relationship between this property and the deletion of classes property of Poisson-Kingman 
models introduced in Pitman (2003). 

Definition 11. [Deletion of classes, Pitman (2003)] Given a random partition II of N, the 
operator deletion of the first k classes is as follows: First let n£ be the restriction of II to 
Hk := N — Gi — • • • — Gk where G\, . . . ,Gk are the first k classes of II in order of their least 
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elements, then derive IL- on N from II£ on H k by renumbering the points of H k in increasing order. 

Proposition 12. The conditional EPPF in Definition 2. of Lijoi, et al. (2008) yields the EPPF of 
a random partition obtained by the operation of deletion of classes. 

Proof. From our group sequential construction the EPPF of 11^ is given by 

/ , v p(m, sin) . , 

psm,n = n ' 1 ' 39 
p{m\n) 

Recall from equation (20) that s < m customers sit at k* new tables in configuration (s±, . . . , Sfc») 
and the remaining m — s customers sit at the old tables in configuration (mi, . . . , mk) for Y2j=i m j = 
m — s, 1 < s < m, y^,j—i Sj — s, rrij > 0, Sj > 1 with probability 

/ \ , I \ p(m + mi > ...,Tik + mk,s 1 ,...,Sk*) ,. n . 

Pm,s(n) = p(m, s n) = , (40) 

p(ni, . . . ,n k ) 



and for Gibbs partitions of type a this yields (cfr. eq. 

k k* 

II^i-aWlICl-a^-iT- (41) 



j'=i i=i 

The denominator in (|39[) is obtained by marginalizing (|40p with respect to s = (si, . . . , Sfc*), i.e. by 
summing over all the ways to partition s observations in k* new tables for every k* 1 i.e. 

I i \ TT/ \ Vn+m,k+k* a — l, a 

p(m|n) = |J_K- -a) m , T r? S Sjk >, , 

i=i fc«=i K "' fc 

hence the EPPF of 11^ is given by 

T/ J« 

/ I \ Vn+m,k+k* I I , , \ 

P( s l m ' n ) = 7} IK 1 " 

Z^fc*=l 'n+ml+t'^s,):* j=l 

which agrees with the result in Proposition 3. in Lijoi et al. (2008). 

Remark 13. Notice that in Lijoi et al. (2008) the result is obtained conditioning on the number 
s of customers in new tables. Given the size m of the new group, this is equivalent to conditioning 
to the number m — s of the observations in old blocks, i.e. to the vector (mi, . . . , rrik) as in our result. 

Now from Proposition 7. in Pitman (2003) if II is a Poisson Kingman PK(p, 7) partition of N, and 
is derived from II by deletion of the first k classes, then IL is a PK(p,j k ) of N where j k is 
given by 7^ = r yQ k , where Q is the Markov transition operator on (0, 00) 

Q(t, dv) = p(t - v)(t - v)t~ 1 f (v)l{0 < v < t}dv 

and f{t) is the probability density corresponding the Levy density p. Since from Gnedin and Pit- 
man (2006) (cfr. item (iii) of Th. 2), we know that an exchangeable partition belongs to the class 
of Gibbs form of type a G (0, 1) if and only if is a mixture of Poisson-Kingman models for some 
mixing density 7, and since IL obtained by deletion of classes of PK(p a ,j) produces PK(p a ,jk) 
which is still of Gibbs form, it should follows that the reproducibility of the Gibbs class also holds 
for infinite conditional structures. This point, which seems to contradict Corollary 2. in Lijoi et al. 
(2008), deserves futher investigation that we postpone to a future paper. 
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